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We study systems with a crossover parameter A, such as the temperature T, which has a threshold 
value A* across which the correlation function changes from exhibiting fixed wavelength (or time 
period) modulations to continuously varying modulation lengths (or times). We report on a new 
exponent, vl, characterizing the universal nature of this crossover. These exponents, similar to 
standard correlation length exponents, are obtained from motion of the poles of the momentum 
(or frequency) space correlation functions in the complex fc-plane (or oj-plane) as the parameter 
A is varied. Near the crossover (i.e., for A — ^ A*), the characteristic modulation wave-vector Ka 
on the variable modulation length "phase" is related to that on the fixed modulation length side, 
q via \Kii — q\ oc \T — T*!"^. We find, in general, that vl = 1/2. In some special instances, vl 
may attain other rational values. We extend this result to general problems in which the eigenvalue 
of an operator or a pole characterizing general response functions may attain a constant real (or 
imaginary) part beyond a particular threshold value. A,. We discuss extensions of this result to 
multiple other arenas. These include the axial next nearest neighbor Ising (ANNNI) model. By 
extending our considerations, we comment on relations pertaining not only to the modulation lengths 
(or times) but also to the standard correlation lengths (or times). We introduce the notion of a 
Josephson timescale. We comment on the presence of aperiodic "chaotic" modulations in "soft-spin" 
and other systems. These relate to glass type features. We discuss applications to Fermi systems - 
with particular application to metal to band insulator transitions, change of Fermi surface topology, 
divergent effective masses, Dirac systems, and topological insulators. Both regular periodic and 
glassy (and spatially chaotic behavior) may be found in strongly correlated electronic systems. 

PACS numbers: 05.50.-|-q, 75.10.Hk, 75.60.Ch 



I. INTRODUCTION 

In complex systems, there are, in general, possibly 
many important length and time scales that characterize 
correlations. Aside from correlation lengths describing 
the exponential decay of correlations, in some materials 
there are length scales that characterize periodic spatial 
modulations or other spatially non-uniform properties as 
in Fig. [TJ We investigate the evolution of these length 
scales as a function of some parameter A. This param- 
eter may be the temperature, the chemical potential, or 
some other physical quantity relevant for description of 
the system being studied. To illustrate our basic premise, 
we will largely focus on temperature dependences of the 
correlation function in this work. However, with a trivial 
change of variables, our results are valid for any parame- 
ter that, when tuned, connects a phase with continuously 
varying modulation lengths (or times) to one in which the 
modulation length (or time) is pinned to a fixed value. 
The crossovers we consider are not symmetry breaking 
transitions. Consequences of our considerations also re- 
late to correlation lengths as we will comment on later. 

Many systems exhibit subtle changes in their corre- 



lation functions at certain special temperatures. We 
focus on situations wherein as the temperature is var- 
ied across a certain crossover temperature T*, an un- 
modulated phase of a system (appearing at T > T.^ 
or at T < T*) may start exhibiting modulations (at 
T < or at T > T* respectively) sans a thermody- 
namic phase transition at T*. A generalization of this 
occurs when modulations in a system are characterized 
by a fixed wavelength on one side of a crossover tem- 
perature and by continuously varying wavelengths on 
the other side. Such an occurrence may generally rear 
its head when interactions of different scales compete 
with one another. A wealth of interesting periodic spa- 
tial patters appear arenas: e.g., the manganites,0 pnic- 
tide d, 0] and cuprate [sl-IToj superconductors, quantum 
Hall systems, [0t3| dense nuclear matter. [l^. [l5[ mag- 
netic systemsjl6|-|2l| heavy fermion compounds, [22|, 
membranes , ^cholesterols , [25| magnetic garnets , 
dipolar systems. (27l[28j systems with nematic phases, 
and countless other systems, jsol - ls^ 



II. OUR MAIN RESULTS AND THEIR 
IMPLICATIONS 
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In this work, we report on the temperature (or other 
parameter) dependence of emergent modulation lengths 
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FIG. 1: Sub-unit-cell resolution image of the electronic struc- 
ture of a cuprate superconductor at the pseudo-gap energy. 
Inset shows Fourier space image of the same figure. Nematic 
and smectic phases are highlighted using the red and blue 
circles respectively. The nematic phase is characterized by 
commensurate wave-vectors Q. The smectic wave-vector, on 
the other hand takes incommensurate values, S which is de- 
endent on the amount of doping, albeit weakly. (From Ref. 
Reprinted with permission from AAAS.) 



that govern the size of various domains present in some 
systems. In its simplest incarnation, our central result 
is that if fixed wavelength modulations characterized by 
a particular finite length scale, L*, appear beyond some 
temperature then, the modulation length, Ld on the 
other side of the crossover differs from as 



\L 



D 



LA (x\T — T^, 



When there are no modulations on one side of 
— )■ oo, we have near the crossover, 



(1) 



(2) 



Apart from some special situations, we find that irre- 
spective of the interaction, = 1/2. We arrive at this 
rather universal result assuming that there is no phase 
transition at the crossover temperature T*. Our result 
holds everywhere inside a given thermodynamic phase of 
a system. 

Our considerations are not limited to continuous 
crossovers. A corollary of our analysis pertains to sys- 
tems with discontinuous ( ''first- order" like) jumps in the 
correlation or modulation lengths. 



We will further comment on situations in which 
wherein a branch point appears at T*. We will present 
examples where we obtain rational and irrational expo- 
nents and also the anomalous critical exponent rj. Our 
analysis affords general connections to the critical scaling 
of correlation lengths in critical phenomena. 

Our results for length scales can be extended to 
timescales. We will, amongst other notions, in employ- 
ing a formal interchange of spatial with temporal coordi- 
nates, introduce the concept of a Josephson timescale. 

Lastly, further deepening the analogy between results 
in the spatial and time domain, we will comment on 
the presence of phases with aperiodic spatial "chaotic" 
modulations (characteristic of amorphous configurations) 
in systems governed by non-linear Euler-Lagrange equa- 
tions. Aperiodic "chaotic" modulations may appear in 
strongly correlated electronic systems. 

In the appendix, we present applications to Fermi 
systems pertaining to metal-band insulator transition, 
change of Fermi surface topology, divergence of effective 
masses, Dirac systems and topological insulators. 



III. THE SYSTEMS OF STUDY 

In this work, we will predominantly consider transla- 
tionally invariant systems on a lattice whose Hamiltonian 
is given by 



H^lY.^i\x~y\)Six)Siy). 



(3) 



The quantities {S'(a?)} portray classical scalar spins or 
fields. The sites x and y lie on a d-dimcnsional hyper- 
cubic (or some other) lattice with N sites. We will set 
the lattice constant to unity. [In the quantum arena, we 
replace the spins S{x) in Eq. by Fermi or Bose or 
quantum spin operators.] 

The results that will be derived in this work apply to a 
variety of systems. These include theories with trivial n- 
component generalizations of Eq. ([3]). In the bulk of this 
work, the Hamiltonian has a bilinear form in the spins. 
We will however, later on, study "soft" spin model with 
explicit finite quartic terms as we now expand on. An 
n-componcnt generalization of Eq. ([3|) is given by the 
Hamiltonian 

H = lY.^{\S-y\)Six)-Siy) + 



-^Y.[s{S)-S{x) 



(4) 



Such a Hamiltonian represents standard (or "hard" ) spin 
or 0{n) systems if m ^ 1 in the large u limit, the quartic 
term enforces a "hard" normalization constraint of the 
particular form S{x) ■ S{x) = n. For finite (or small) m. 
Equation (j?]) describes "soft" -spin systems wherein the 
normalization constraint is not strictly enforced. 
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In what follows, i;(fc) and s{k) will denote the Fourier 
transforms of V^(|x — j7|) and S{x). We employ the fol- 
lowing Fourier conventions, 



a{k) = ^A(f)e'* 

X 



— ik-: 



(5) 



With these conventions in tow, in Fourier space, Eq. ([3]) 
reads 



(6) 



When v{k) is analytic in all momentum space coordi- 
nates, it is a function of |fcp = (and not a general 

function of k ; 



^^i^i kf with {ki} being the Cartesian 

components of k). This is so as |fc| has branch cuts when 
viewed as a function of a particular fc; (with all other 
ki'^i held fixed). The lattice Laplacian that links nearest 
neighbors sites in real space becomes 



d 

(E 

1=1 



2(V(l-cosfcO 



(7) 



in fc-space. A^: veers towards |/cp in the continuum (small 
k) limit. The two point correlation function for the sys- 
tem in Eq. ^ is, G{x) = {S{0)S(x)). At large distances 
X = the correlation function has a general asymptotic 
behavior 




(8) 



In the i-th term, fi{x) is an algebraic prefactor, is 
the modulation length and is the corresponding corre- 
lation length. In general, the function fi{x) may contain 
a factor with an anomalous exponent (usually not an 
integer), such as, fi{x) oc l/a;'^^^+''. Generally, there 
can be multiple correlation and modulation lengths. In 
Fourier space, Gik) = ^(|s(fc)|2). The modulation and 
correlation lengths can be obtained respectively from the 
real and imaginary parts of the poles of G(fc) in the com- 
plex fc-plane. 



General considerations: Correlation and modulation 
lengths from momentum space correlation function 

The correlation function G{x) in (d-dimcnsional) real 
space is related to the momentum space correlation func- 
tion G{k) by 



G{x) 



df^k 



G{k)e 



-ik-x 



(9) 



On the lattice, the integral above must be replaced by 
summation over fc- values belonging to the first Brillouin 
zone. In the continuum, which we discuss here, the inte- 
gral range is unbounded. Even in lattice systems, doing 
an unbounded summation over fc- values provides a good 
approximation for the correlation function in real space 
in many scenarios. 

For spherically symmetric problems, i.e., when G{k) ~ 
G{k), 



G{x) 



k 



dfcJ./.-i(M^(^^^ (10) 



(27r)'^/2 (fca;)^/2 



where J,y(x) is a Bessel function of order v. The above 
integral can be evaluated by choosing an appropriate con- 
tour in the complex fc-plane. The contour can be closed 
along a circular arc of radius i? — ?> oo provided 



|G(fc)|<fc- 



0jS k 



(11) 



In evaluating the integral in Eq. (|10p . we obtain con- 
tributions from residues associated with the poles of the 
integrand as well as contributions from its branch points. 
We use K = + iKj to represent the poles and branch 
points of the integrand in the complex plane. The cor- 
relation and modulation lengths in the system are deter- 
mined respectively by the imaginary {Kj) and real parts 
(Kr) of these poles and branch points. Together, all 
these singularities can be compactly expressed as 



1 



G(™)(A') 



(12) 



where < m < oo is the order of the smallest order 
derivative of G'(fc) which diverges at fc = i^.jlHl 

In footnote [36| , we comment on the situation in which 
the function G(T, fc) is an entire function of fc (i.e., when 
is G is analytic everywhere). 



IV. A UNIVERSAL DOMAIN LENGTH 
EXPONENT - DETAILS OF ANALYSIS 

We now derive (via various inter- related approaches), 
our central result - the existence of a new exponent for 
the domain length in rather general systems with real or 
complex scalar fields, vectorial (or tensorial) fields of both 
the discrete (e.g.. Potts like) and continuous variants. 

We will now consider the situation in which the sys- 
tem exhibits modulations at a fixed wave-vector q for a 
finite range of temperatures on one side of T*, [viz., (i) 
T > Ti,, or, (ii) T < T*] and starts to exhibit variable 
wavelength modulations on the other side [(iii) T < T^, 
for (i) and T > T^, for (ii)]. A schematic illustrating 
this is shown in Fig. [2l In sub-section IIV Al we will as- 
sume that the pair correlation function is meromorphic 
(realized physically by absence of phase transitions) at 
the crossover point and illustrate how modulation length 
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(a) Variable modulation length. (b)The crossover point (T = T, ). (c)Fixcd modulation length. 



FIG. 2: Schematic showing the trajectories of the singularities of the correlation function near a fixed - variable modulation 
length crossover. Two poles of the correlation function merge at = fc* at T = T*. On the fixed modulation length side of the 
crossover point, Re k = q. 



exponents may appear. In sub-section IIV Bl we will com- 
ment on the situation where the crossover point may be 
a branch point of the correlation function. 



A. Crossovers at general points in the complex 
fc-plane 

In the up and coming, we will assume that the pair 
correlator, G{T, k) is a meromorphic function of k and 
T near a crossover point. Our analysis below is exact 
as long as we do not cross any phase boundaries. Such 
a case is indeed materialized in the incommensurate- 
commensurate crossovers in the three-dimensional ax- 
ial next-nearest-neighbor Ising (ANNNI) model [13, l38l | 
(which is of type (ii) in the classification above). This 
phenomenon is also seen in the ground state phase dia- 
gram of Frenkel-Kontorova models [s^ in which one of 
the coupling constants is tuned instead of temperature. 

In the following, we present two alternative derivations 
for the universal exponent characterizing this crossover. 



1. First approach 

In general, if the pair correlation function G{T, k) 
is a meromorphic function of the temperature T and 
the wave- vector k near a crossover point (T^fc*), then 
G~^{T, k) must have a Taylor series expansion about that 
point. We have. 



G-i(T,fc) 



OO 

E 



Ar. 



,(r-r,)'"i(fc-fc,)™^(i3) 



mi ,7712—0 



Since G~^(r*,fc*) = 0, we have, Aqo = 0. Let us try to 
find the trajectory of the pole K{T) (with KiT^,) = fc*) 
of G(T, k) in the complex fc-plane as the temperature is 
varied around T*. Writing down the leading terms of 



G ^{T, k) , we have, in general, 

[m/a] 

G-^T, fc) - ^{T- T,)™-'^P(fc - + 

o((r-T,)™(fc-fc,)"), (14) 

as (T, k) — (T*, kif) with m, n, a, b integers, m, n > and 
a,b > 1, [x] represents the greatest integer less than or 
equal to X and o{x) represents terms negligible compared 
to X. We have, 



K{T) - k,+C{T~T,) 



a/b 



(15) 



where C is some constant, yielding vl = a/b. By the 
very definition of T*, on one side of T* [(i) or (ii) above], 
there exists at least one root [say, K{T)] of G~^ satisfying 
Kr{T) = q, where g is a constant. On the other side [(iii) 
above], Kpi{T) ^ q. As such, the function K{T) is non- 
analytic at . The left hand side of Eq. ([T5|) is therefore 
not analytic at T = , implying that the right hand side 
cannot be analytic. This means that (a/b) cannot be an 
integer, which in turn implies that b > 2. Therefore, in 
the most common situations we might encounter, 



G-^{T,k) - A{T- 
a = \ and 6 = 2. 



T^) + B{k~Kf 



(16) 



When Fourier transforming G{T, k) by evaluating the in- 
tegral in Eqs. ([ni[TU]) using the technique of residues, the 
real part of the poles (i.e., Kpt) gives rise to oscillatory 
modulations of length Lo = Kn- If the modulation 
length locks its value to 27r/g on one side of the crossover 
point, then, on the other side, near T*, it must behave as 



\2n/LD-q\ oc |T-T*|i/2 



= 1/2. 



(17) 



2. Second approach 



We now turn to a related alternative approach that 
similarly highlights the universal character of the modu- 
lation length exponent. If the correlation function G is 
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FIG. 3: Location of the poles of the correlation function of 
the large n Coulomb frustrated ferromagnet for J — Q — 1 
in the complex fc-plane. The circle and the K-axis show the 
trajectory K(T) of the poles as the temperature T is varied. 



a meromorphic function of k, then, expanding about a 
zero Ki{T) of G^^ , we have, 

G-\T, k) = A{T) [k - K,{T)r' G^\T, k), (18) 

where Gi^{T,k) is an analytic function of k and 
G];'^{T,Ki{T)) ^ 0. We can do this again for the func- 
tion Gi^{T, k) choosing one of its zeros K2{T) and con- 
tinue the process until the function left over does not 
have any more zeros. We have, 

G-i(r, k) = AiT) n (fc - i^a(T))™° G;\T, fc), (19) 

a=l 

where the function G~^{T, k) is an analytic function with 
no zeros, niaS are integers and, in principle, p may be 
arbitrarily high. This factorization can be done in each 
phase where G is meromorphic. Let Ki{T) be a non- 
analytic zero of G~^, i.e., one for which Re Ki{T) = q 
on one side of T = T.^ . To ensure analyticity of G~ ^ in 
T in the vicinity of T = T*, there must be at least one 
other root K2{T), such that as T T*, both Ki{T) and 
K2{T) veer towards fc*, where Re fc* = q [e.g., see Fig. 
[3] which is of type (i) above, fc* ~ ±i]. In other words, 
p in Eq. cannot be smaller than two. The proof of 
this assertion is simple, li p = 1, then, according to Eq. 
(Ull), G-\T,k) = A{k - R\{T))G^\T,k). At T = 
however, Ki{T) is not analytic, implying that G~^{T, k) 
can be analytic only if p > 2. For p > 2, at T*, 



will, to leading order, vary quadratically in (fc — fc*) in 
the complex k plane near fc*. Thus, 



dG- 



dk 



(20) 



{T,,k,) 



Now, if G~^ has a finite first partial derivative relative to 
the temperature T then, for a pole K near fc* , to leading 
order. 



G-^(r*,fc*) + (T-T.,) 



dG- 



dT 

{K-k,f d^G-' 



2! 



{T,,k.) 



0. (21) 



By its definition, fc, satisfies the equality G ^ (T* , k^ ) 
0. Therefore, 



\K-kJ = 



\ 



2{n - T) 



dG~ 



dT 



{T„k,) 









(T,.k,) 



(22) 



Equation (|17p is an exact equality. It demonstrates 
that the exponent ~ 1/2 universally unless one 



and vanishes at (T*,fc*).[4^ Often, 



G ^(T, k) is a rational function of fc, i.e. 



G-\T,k) 



G-HT,k) 

G-\T,ky 



(23) 



where G^^{T, k) and G^^{T, k) are polynomial functions 
of fc. In those instances, we get the same result as above 
by using G~^{T,k) in the above arguments. The value 
of the critical exponent is similar to that appearing for 
the correlation length exponent for mean-field or large n 
theories. It should be stressed that our result of Eq. P7)) 
is far more general. 

Lock-in of the correlation length. Apart from the 
crossovers across which the modulation length locks in to 
a fixed value, we can also have situations where the cor- 
relation length becomes constant as a crossover temper- 
ature T** is crossed. If this happens, our earlier analysis 
for the modulation length may be replicated anew for the 
correlation length. Therefore, if the correlation length 
has a fixed value on one side (T < T** or T > T**) of 
the crossover point, then, on the other side (T > T** or 
T < T**, respectively), near T**, it must behave as. 



|i/e-i/Co|(x|r-r„r% 



(24) 



where, like vl, = 1/2 apart from special situations 
where it may take some other rational values. Here and 
throughout, we use v,. to represent the usual correlation 
length exponent, v to distinguish it from the modulation 
length exponent vj^. 



6 



B. Branch points 

A general treatment of a situation in which the 
crossover point is a branch point of the inverse correla- 
tion function in the complex fc-plane is beyond the scope 
of this work. Branch points arc ubiquitous in correlation 
functions in both classical as well as quantum systems. 

For example, in the large n rendition of a bosonic sys- 
tem (with a Hamiltonian of Eq. ^ and S{x) repre- 
senting bosonic fields), the momentum space correlation 
function at temperature T is given by [ill, |42| 



G{k) 



v{k) + 



JIB 



,(25) 



where fii is a constant having dimensions of energy, ^ is 
the chemical potential, nsix) = l/(e^ — 1) is the Bose 
distribution function and fc^ is the Boltzmann's constant. 

Similar forms, also including spatial modulations in 
G{r), may also appear. We briefly discuss examples 
where we have a branch cut in the complex /c-plane. 

The one-dimensional momentum space correlation 
function. 



G{k) 



1 



1 



^{k - qY + r ^{k + qY+r 
reflects a real space correlation function given by 



(26) 



G{x) 



2 cos((7x)Ko(x-y/r) 



where Ko(-) is a modified Bessel function. Thus, as is to 
be expected, we obtain length scales associated with the 
branch points K = ±q ± i^/r. 

Similarly, the three-dimensional real space correlation 
function corresponding to 



G{k) 



1 



^{k-qf+r 



(28) 



exhibits the same correlation and modulation lengths 
along with an algebraically decaying term for large sepa- 
rations. Another related G~^(fc) involving a function of 

|fc| (i.e., not an analytic function of fc^) was investigated 
earlier. [4J] 

Throughout the bulk of our work, we consider simple 
exponents associated with analytic crossovers. In con- 
sidering brach points, our analysis may be extended to 
critical points. As is well known, at critical points of d 
dimensional systems, the correlation function for large r, 
scales as 



G{r) cx 



1 



(29) 



with the anomalous exponent. Such a scaling implies, 
for non-integer 77, the existence of a A: = branch point 
of G(A:). 



If the leading order behavior of 1/ G^™^ (T, k) is alge- 
braic near a branch point (T* , ^* ) , then we get an alge- 
braic exponent characterizing a crossover at this point [m 
being the lowest order derivative of G(fc) which diverges 
at fc = fc* as in Eq. (|T2|) ]. That is, we have, 

5^ A{T -T,y^ - B{k-KY^ 



G('^){T, k) 



as (T,fc) ^ (T,,fc*). 



(30) 



This implies that the branch points K deviate from fc* 
as 

[K-K)^^-^ {T-nr'^\ (31) 

We therefore observe a length scale exponent v = zxjzi 
at this crossover. This exponent may characterize a cor- 
relation lengths and/or a modulation lengths. The ex- 
ponent z\l z-2, may assume irrational values in many sit- 
uations in which the function G~^(T,k) is not analytic 
near the crossover point. Such a situation could give rise 
to phenomena exhibiting anomalous exponents 77. For 
example, if we have a diverging correlation length at a 
critical temperature Tc, for a system with a correlation 
function which behaves as in Eq. (|29p . then, we have in 
Eq. dSni), Z2 = 2 - 77. Thus, we have. 



Ld - Ldc\ (x\T -Tc 



zi 



(32) 



(27) where Ldc 



27r/|Rc fc*|, and more importantly. 



Zl 



(33) 



Other critical exponents could also, in principle, be cal- 
culated using hyper-scaling relations. 

If G~^{T,k) has a Puiseux representation about the 
crossover point, i.e.. 



G-i(r,fc) 



00 00 

E 

m=mo p=Po 



amp{k ~ k,r^-{T - T,)P/^(34) 



with Omopo = 0, where mg, po, a and b are integers, 
then, the result we derived above applies to the relevant 
length scale and the crossover exponent v = a/b^ is again 
a rational number. 

Generalizing, if G^^{T, k) is the ratio of two Puiseux 
series, we use the numerator to obtain the leading order 
asymptotic behavior and hence obtain a rational expo- 
nent. 



C. A corollary: Discontinuity in modulation 
lengths implies a thermodynamic phase transition 

Non-analyticities in the correlation function G(fc) for 
real wave- vector k imply the existence of a phase transi- 
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tion. This leads to simple corollaries as we now briefly 
elaborate on. A sharp discontinuous jump in the value 
of the modulation lengths (and/or correlation lengths) 
implies that the zeros {Ka} of G~^{k) in the complex 
k plane, exhibit discontinuous ("first order-like") jumps 
as a function of some parameter (such as the tempera- 
ture T when T = T*). When this occurs, as seen by, 
e.g., differentiating the reciprocal of the product of Eq. 
(|T9|) . the correlation function will, generally, not be an- 
alytic as a function of T at T = T*. Putting all of the 
pieces together, we see that a discontinuous change in the 
modulation (or correlation) lengths impies the existence 
of a bona fide phase transition. Thus, all commensurate- 
commensurate crossovers must correspond to phase tran- 
sitions. For example, see the ANNNI model. |44l| 

D. Diverging correlation length at a spinodal 

transition 

Our analysis is valid for both annealed and quenched 
systems so long as translational symmetry is main- 
tained (and thus, the correlation function is diagonal 
in fc-space). In particular, whenever phase transitions 
are "avoided" the rational exponents of Eq. ((T5|) will 
appear.gllilliil 

In diverse arenas, we may come across situations in 
which there are no diverging correlation lengths even 
when the inverse correlation function has zeros corre- 
sponding to real values of the wave- vector. These are 
signatures of a first order phase transition, e.g., transi- 
tion from a liquid to a crystal. If the first order phase 
transition is somehow avoided, then the system may enter 
a metastable phase and may further reach a point where 
the correlation length diverges, e.g., a spinodal point. If 
it is possible to reach this point and if the inverse cor- 
relation function is analytic there, then our analysis will 
be valid, thereby leading to rational exponents charac- 
terizing the divergence of the correlation length. There 
are existing works in the literature which seem to suggest 
that such a point may not be reachable. For example, in 
mode coupling theories of the glass transition, the sys- 
tem reaches the mode coupling transition temperature 
Tmct at which the viscosity and relaxation times diverge 
and hence does not reach the point where the correlation 
length blows up.f47| 

E. Conservation of characteristic length scales 

In Ref. it was mentioned that the total num- 

ber of characteristic length scales in a large-n system re- 
mains constant in systems in which the Fourier space 
interaction kernel v{k) is a rational function of and 
the real space kernel is rotationally invariant. (Simi- 
lar results hold for systems with reflection point group 
symmetry. [4^ ) In this sub-section, we generalize that ar- 
gument and say that whenever the Fourier space corre- 



lator G{k) of a general rotationally invariant system is a 
rational function of fc^,i.e., 

the total number of correlation and modulation lengths 
remains constant apart from isolated points as a tuning 
parameter A is smoothly varied. In Eq. 1351 the func- 
tions P(fc^) and Q{k^) are polynomial functions of k^. 
Rotational invariance requires that G{k) is real-valued 
for real wavevectors k. As argued in Ref. [l^, all length 
scales in the such systems are associated with the poles of 
G{k) in the complex fc-plane and these remain constant 
for a given form of the function G{k). Each real root of 
the function Q{k^) gives rise to a term in the real space 
correlation function which has one correlation or modu- 
lation length. Non-real roots (which necessarily come in 
complex conjugate pairs) give rise to a correlation and a 
modulation length. Thus, the total number of character- 
istic length scales in the system is equal to the order of 
the polynomial function Q{k^) which remains fixed. 

V. 0(n) SYSTEMS 

The correlation function for 0{n) systems can be cal- 
culated exactly at both the low and the high temperature 
limits. At intermediate temperatures, various crossovers 
and phase transitions may appear. In this section, we 
discuss the low and high temperature behavior length 
scales characterizing 0{n) systems. 

A. Low temperature configurations 

It was earlier demonstrated [i^ that for 0{n > 2), all 
ground states of a system have to be spirals (or poly- 
spirals) of characteristic wave- vectors Qa, given by 

v{qa) — — min v{k), (36) 

where M'' represents the set of all d-dimcnsional real vec- 
tors. At r = 0, the modulation lengths in the system are 
given by 

L'^"(r = 0) = 27rM,„, (37) 

where i{l < i < d) labels the Cartesian directions in 
d dimensions. This, together with Eq. (|38p gives us 
the high and low temperature forms of the correlation 
fimction and its associated length scales. 

B. High temperatures 

As is well appreciated, diverse s;Y^tems behave in the 
same way at high temperatures. [50[ For 0{n) systems 
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[53 (any n), 

G-\T,k) ^l + v{k)/kBT + 0{l/T'^). (38) 

The high temperature series may be extended and ap- 
phcd at the crossover temperature T* , if there is no phase 
transition at temperatures above and for aU relevant 
real fc's, \v{k) \ ^ ksT^,. [A detailed example will be stud- 
ied in Sec. IV El ] Generally, Eq. ([38]) may be analytically 
continued for complex fc's and in the vicinity of T*, 



Sk 



i\ kBjT.-T) 
w(™)(/c*) 



(39) 



where k^ is a characteristic wave- vector at Tj^. In the 
above, Sk denotes the change in the location of the poles 
K of G^^ when the temperature is changed from to 
T (i.e., 6k = K — k^,) and m is the order of the lowest 
non- vanishing derivative of v{k) at fc*. As in previous 
analysis, v'{k^) = and m > 2. For general u(fc), typi- 
cally m = 2 and i^l = 1/2 as before. 

We now turn to examples which explicitly illustrate 
how our results are realized including exceptional systems 
with non-trivial exponents. 



C. Large n Coulomb frustrated ferromagnet — 
modulation length exponent at the crossover 
temperature T, 

In the current sub-section and the two that follow, we 
will discuss the large n limit in 0{n) systems. The re- 
sults in the previous two sections pertain to arbitrary n. 
We illustrate how our result applies to the large n 1511 
Coulomb frustrated ferromagnet. As is well known j52| . 
in the large n limit, 0{n) systems are exactly solvable 
and behaves as the spherical model, [s^ The correlation 
function in fc-space is given by 



G~\T,k) = [v{k)+fi{T)]/kBT, 



(40) 



where v{k) is the Fourier space interaction kernel and 
^(T) is a Lagrange multiplier, see e.g. Ref. [iil,!!^], that 
enforces the spherical constraint 



of strength Q) is given by v{k) = Jk"^ + Q/k"^, where 
J and Q are positive constants. The critical tempera- 
ture, Tc of this system is given by ii{Tc) = —2^JQ. 
At Tcj the correlation length is infinity and the modu- 
lation length IS Ld = 2tt ^ J/Q. As the temperature is 
increased, the modulation length increases and the corre- 
lation length decreases. At T*, given by ^{T^,) = 2y/JQ^ 
the modulation length diverges and the correlation length 
becomes ^ = {/ J IQ- At temperatures above T, , the cor- 
relation function exhibits no modulations and there is one 
decreasing correlation length and one increasing correla- 
tion length. The term in the correlation function with 
the increasing correlation length becomes irrelevant at 
high temperatures because of an algebraically decaying 
prefactor. The divergence of the modulation length at T* 
shows an exponent of z^l = 1/2. ji^ 



D. An example with 7^ 1/2 

In what follows, we demonstrate, as a matter of prin- 
ciple, that the exponent for the divergence of the mod- 
ulation length (and also the correlation length) can be 
different from 1/2 in certain special cases. As an illustra- 
tive example, we consider a large n (or spherical model) 
system for which in Eq. 

v{k) = ^(fc2 + /72^2 ^4S(fc2 + ^-2) 

^^CIOe^l-^) + Dl(k^ + l-^)\ (42) 

where Is is a screening length. If we set A = B = 
C = D ^ 1 then in the resultant system jy^, 7^ 1/2 at a 
crossover temperature. It has a critical temperature Tc, 
given by fJ.{Tc) = —10. At Tc, the modulation length is 



l^(5(f).s'(f)) = l. 



(41) 



Ld = 27r/-\/l — 1/^2 and the correlation length blows up 
(as required by definition). At the crossover temperature 
(for which /i(r*) = 6) the modulation length diverges 
and the correlation length scales as ^ = 1/ -^/l + l/If A 
temperatures just below T*, the modulation length Ld 
diverges as Ljj oc {T^ — T)"^/"* meaning that i/^ = 1/4. 
This is because the first three derivatives of v{k) vanish 
at fc = i, which is the characteristic wave- vector at 
(see Fig. 15). 



E. An example in which T, is a high temperature 



The paramagnetic transition temperature Tc is obtained 
from the relation, fJ-{Tc) = — min^gR ^(fc). Below Tc, 
the Lagrange multiplier fi{T) = jiiTc)- Above Tc, 
fJ-{T) is determined by the global average constraint that 
G{x = 0) = J2k^(^) = 1- This global constraint also 
implies that, above Tc, small changes in temperature re- 
sult in proportional changes in fi{T) and at high tempera- 
tures, /x(T) is a monotonic increasing function of T. The 
Fourier space kernel v{k) for the "Coulomb frustrated 
ferromagnet" (in which nearest neighbor ferromagnetic 
interactions of strength J compete with Coulomb effects 



We now provide an example in which the high tem- 
perature result of Sec. IV Bl (valid for any 0(n) system 
with arbitrary n) can be applied at a crossover point. 
Consider the large n system in Eq. ([42]) with A = 1, 
B ^ 1, C — 1/4, D = and the screening length. 
Is ^ B. The critical temperature of this system is 
given by fJ,{Tc) ^ — 4-\/B where the modulation length 
is Ld ^ 27r-v^4B. There is a crossover temperature 
such that ^{Tif) ^ AB^ . One of the modulation lengths 
diverges at T^. The corresponding correlation length is 
given by ^ Xj \J2B. This provides an example in which 
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FIG. 4: Location of the poles of the correlation function of 
the system in Eq. (|42|l for large Is (small screening) in the 
complex fc-plane. 



\v{k) \ <^ ksT^ for all real fc's satisfying < tt. The sec- 
ond derivative of v{k) is non-zero at the crossover point, 
resulting in a crossover exponent = 1/2. 



VI. CROSSOVERS IN THE ANNNI MODEL 

We now comment on one of the oldest studied examples 
of a system with a crossover temperature. The following 
Hamiltonian represents the ANNNI model. jSTllaslliil 



H = -,hY.S{x)S{y) + J2 S{x)S{y), 



(43) 



where as throughout, x is a lattice site on a cubic lattice, 
and the spins S(x) ~ ±1. The couplings, Ji, J2 > 0. In 
the summand, (•) represents nearest neighbors and ((•)) 
represents next nearest neighbors along one axis (say the 
Z-axis), see Fig. [S] Depending on the relative strengths 
of Ji and J2, the ground state may be either ferromag- 
netic or in the "(2) phase". The "(2) phase" is a pe- 
riodic layered phase, in which there are layers of width 
two lattice constants of 'up" spins alternating with lay- 
ers of "down" spins of the same width, along the Z-axis. 
As the temperature is increased, the correlation function 
exhibits jumps in the modulation wave-vector at differ- 
ent temperatures. At these temperatures, the system 
undergoes first order transitions to different commensu- 
rate phases. The inverse correlation function G~^{T,k) 




FIG. 5: The coupling constants in the three-dimensional 
ANNNI model. 
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FIG. 6: Existence of incommensurate phases between the 
commensurate regions in the phase diagram of the ANNNI 
model, (a) Mean field phase diagram of the ANNNI model in 



three dimensions. The shaded regions show higher order com- 
mensurate phases which have variable modulation length in- 
commensurate phases in between (From Ref. [sj. Reprinted 
with permission from APS.) (b) Phase diagram for the three- 



dimensional ANNNI model using a modified tensor product 
variational approach (From Ref. [s^. Reprinted with per- 
mission from APS.) [(c)] Variation of wavelength along paths 
AiBi and A3B3 of |(b)| showing a smooth variation of the 
wavelength near the paramagnetic transition line (From Ref. 
[55l |. Reprinted with permission from APS.) (d) Cartoon of 



an incommensurate-commensurate crossover region from (a) 



is therefore not an analytic function of k and T at the 
transition points. The phase diagram for the ANNNI 
model, however, also has several crossovers where the sys- 
tem goes from a commensurate phase to an incommensu- 
rate phase with a continuously varying modulation length 
(see Fig. inj.jl^, At these crossover points, follow- 
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ing our rigorous analysis, we expect a crossover exponent 
vl = 1/2- Such a scaling of the modulation length has 
been estimated by several approximate techniques near 
the "Lifshitz point" Pl . [M [5i-l6l| The Lifshitz point is 
the point in the phase diagram of the ANNNI model at 
which the high temperature paramagnetic phase coex- 
ists with the ferromagnetic phase as well as a phase with 
continuously varying modulation lengths. It is marked as 
Pl in Fig. |6(b)[ Although the point Pl has a first order 
transition, it can be thought of as a limit in which the 
incommensurate and commensurate regions in Fig. 6(a) 
shrink and merge to a single point. We would also like 
to point out that it is known |62| that if the wave-vector 
takes all possible rational values ("complete devil's stair- 
case"), we have no first order transitions. Additionally, 
non-analyticity of the correlation function does not pro- 
hibit other quantities from having continuous crossover 
behavior. For example, the correlation of the fluctua- 
tions, i.e., the connected correlation function may gener- 
ally exhibit continuous variation from a fixed to a variable 
modulation length phase. If the inverse connected corre- 
lation function is analytic, our result can be applied to it 
resulting in a crossover exponent of 1/2. 



VII. PARAMETER EXTENSIONS AND 
GENERALIZATIONS 

It is illuminating to consider simple generalizations of 
our result to other arenas. We may also replicate the 
above derivation for a system in which, instead of tem- 
perature, some applied other field A is responsible for the 
changes in the correlation function of the system. Some 
examples could be pressure, applied magnetic field and 
so on. The complex wave- vector k could also be replaced 
by a frequency w whose imaginary part would then cor- 
respond to some decay constant in the time domain. 

More generally, we look for solutions to the equation 



G-i(A,ii) =0, 



(44) 



with the variable u being a variable Cartesian component 
of the wave-vector, the frequency, or any other momen- 
tum space coordinate appearing in the correlation func- 
tion between two fields (u = ki,Lu, and so on). Replicat- 
ing our steps mutatis mutandis, we find that the zeros of 
Eq. scale as |w — uqI oc |A— A*|^/^ whenever the real 
(or imaginary) part of some root becomes constant as A 
crosses A*. Thus, our predicted exponent of vl = 1/2 
could be observed in a vast variety of systems in which a 
crossover occurs as the applied field crosses a particular 
value, in the complex wave-vector like variable. 

Another generalization of our result proceeds as 
follows. [63| Suppose that we have a general analytic op- 
erator (including any inverse propagator) G~^(A) that 
depends on a parameter A. Let be a particular eigen- 
value, 



The secular equation for the eigenvalues of is an an- 
alytic function in A. We may thus replicate our earlier 
considerations to obtain similar results. In doing so, we 
see that if aa{X) changes from being purely real to be- 
coming complex as we vary the parameter A beyond a 
particular threshold value A* (i.e., if aa{X > A*) is real 
and aa{X < A*) is complex, or the vice versa), then the 
imaginary part of aa{X) will scale (for A < A* in the 
first case noted above and for A > A* in the second one) 
as Im {act (A)} cx |A — A*|^/^. A particular such realiza- 
tion is associated with the spectrum of a non-Hermitian 
Hamiltonian [playing the role of in Eq. (05])] which, 
albeit being non-Hermitian, may have real eigenvalues 
(as in VT symmetric Hamiltonians) . [gj] In this case, the 
crossover occurs when a system becomes VT symmetric 
as a parameter A crosses a threshold A* . 

Similarly, if aa{X) changes from being pure imaginary 
to complex at A = A*, then the real part of the eigenvalue 
will scale in the same way. That is, in the latter instance. 
Re {a„(A)} cx |A- A^l^/^. 

Our next brief remark pertains to some theories with 
multi-component fields, e.g^n component theories with 
Hamiltonians of the form. 



41| 



H 



2N ^ 



Vij{k)si{k)sj{k), 



(46) 



k,i,j 



G-i(A)|a„(A)) =a„(A)|a„(A); 



(45) 



in which, unlike Eq. ([6]) (as well as standard 0{n) theo- 
ries), the interaction kernel Vij might not be diagonal in 
the internal field indices i,j = l,2,...,n. An example 
is afforded by a field theory in which n component fields 
are coupled minimally to a spatially uniform (and thus 
translationally invariant) non-Abelian g aug e background 
which emulates a curved space metric. [4 1| In this case, 
the index a in Eq. (|45p is a composite of an internal field 
component coordinate w — 1, 2, . . . , n and A:-space coor- 
dinates. For each of the n branches w, we may determine 
the associated fc-space zero eigenvalue of Eq. (|45p which 
we label by (i.e., au,,k=K^{X) = 0). The largest corre- 
lation is length is associated with the eigenvector which 
exhibits the smallest value of |Im K^\. As usual, as A 
is varied, we may track for each of the n branches, the 
trajectories poles of G in the complex fc-plane. Although 
the location of the multiple poles may vary continuously 
with the parameter A, the dominant poles (those associ- 
ated with the largest correlation length) might discontin- 
uously change from one particular subset of eigenvectors 
to another (see Fig. [7]). As such, the correlation function 
of the system may show jumps in its dominant modu- 
lation length at large distances as A crosses a thresh- 
old value A* even though no transitions (nor cross-overs 
similar to that of Fig. © which form the focus of this 
work) are occurring. Such jumps in the large distance 
modulation lengths appear in 0(n) systems defined on a 
fixed, translationally invariant, non-Abelian background 
or metric as in Ref. [4lj . 

In Appendix [XI we discuss exponents associated with 
lock-ins of correlation and modulation lengths in Fermi 
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FIG. 7: "Jumps" in the modulation length: The figure shows 
the evolution of the poles associated with two different eigen- 
vectors with the parameter A in the complex fc-plane. The 
solid portion of the trajectories show which pole corresponds 
to the dominant term (larger correlation length) in the cor- 
relation function. The x-s denote the poles at A = A» and 
the arrows denote the direction of increasing A. It is evi- 
dent, therefore, that the modulation length corresponding to 
the dominant term jumps from Ldi to Ld2 as A crosses the 
threshold value A, . 



systems. When dealing with zero temperature behavior, 
we use the chemical potential fi as the control parameter 
A. We discuss metal-insulator transition, exponents in 
Dirac systems and topological insulators. Additionally, 
we comment on crossovers related to changes in the Fermi 
surface topology as well as those related to situations 
with divergent effective mass. 



VIII. IMPLICATIONS FOR THE TIME 
DOMAIN: JOSEPHSON TIME SCALES AND 
RESONANCE LIFETIMES 



As we alluded to above, the results that we derived 
earlier that pertained to length scales can also be ap- 
plied to time scales in which case we look at a temporal 
correlation function characterized by decay times (corre- 
sponding to correlation lengths) and oscillation periods 
(corresponding to modulation lengths). We may obtain 
decay time and oscillation period exponents whenever 
one of these time scales freezes to a constant value as 
some parameter A crosses a threshold value A* . 

Many other aspects associated with length scales have 
analogs in the temporal regime. Towards this end, in 
what follows, we advance the notion of a "Josephson time 
scale" . We first very briefly review below the concept of 
a Josephson length scale. In many systems [with correla- 
tion functions similar to Eq. ([2^ ]. just below the critical 
temperature, the correlation function as a function of 



G{k) cx 



for k > I/O, 
for k < I/O, 



(47) 



thus defining the Josephson length scale, 0-|65| Such an 
argument may be extended to a time scale, tj (real or 
imaginary) in systems with Lorentz invariant propaga- 
tors. For a given wave-vector fc, tj may be defined as. 



G{k,Lu) 



-2 



for uj > 1/tj, 
for UJ < 1/tj, 



(48) 



where uj is the frequency conjugate to time while perform- 
ing the Fourier transform and rit{^ 0) is an anomalous 
exponent for the time variable. 

We next briefly allude to another possible simple ap- 
plication of our result. As is well known in high energy 
(see, e.g., Ref. [6^ for a standard textbook treatment) 
and many body theories, the Fourier transform of the two 
two point correlation function G(k,u;) generally exhibits 
isolated poles corresponding to the one particle states 
as well as bound states and a branch cut that reflects 
a continuum of multi-particle states (i.e., two particles 
or more). Such a continuum of states arises when the 
squared four-momentum = j (? — fp' exceeds the 
threshold necessary for the production of two particles, 
i.e., > (2m)^c^ with m the particle rest mass and c 
the speed of light. Single particle (and bound) states and 
continuous multi-particle states lead to the aforemen- 
tioned respective single poles and branch cuts along the 
real axis. We may consider an application of our ideas 
in the vicinity of zero energy bound states (as in, e.g^ the 
Feshbach resonance of the BCS to BEC crossover . j67H70| 
in dilute gases where the crossover is driven by vary- 
ing an attractive contact interaction of strength g) when 
poles on the real axis are just about to splinter into poles 
with a infinitesimal imaginary part. Generally, when, by 
virtue of self-energy corrections, the poles attain a finite 
imaginary part in the p^ plane, the corresponding states 
attain a finite lifetime (with the lifetime being the analog 
of the correlation length/time in the two-point correla- 
tion functions that we discussed hitherto). The relations 
(and exponents) that we derived thus far may be applied, 
mutatis mutandis, for the description of processes asso- 
ciated with the depinning of the poles off the real axis, 
due to the imaginary part of the self energy E, leading 
to resonances with a finite life-time. This relates to the 
scaling of the lifetime r of resonances in cold atomic gases 
as a function of (go — g) where go is the strength of the 
contact interaction at the BCS to BEC crossover point. 



IX. CHAOS AND CLASSINESS 

Thus far, we have considered phases in which the mod- 
ulation length is well defined. For completeness, in this 
section, we mention situations in which aperiodic phases 
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may be observed. The general possibility of such phe- 
nomena in diverse arenas is well known. [62, [ill We focus 
here on translationally invariant systems of the form of 
Eqs. pi4p with competing interactions on different scales 
that lead to kernels such as 



v{k) 



Cik~ 



C2, 



(49) 



where Ci and C2 are positive constants may give rise to 
glassy structures for non zero u. Such a dispersion may 
arise in the continuum (or small k) limit of hyper-cubic 
lattice systems with next nearest neighbor interactions 
(giving rise to the /c^ term) and nearest neighbor in- 
teractions (giving rise to the k'^ term). Within replica 
type approximations, such kernels that have a finite k 
minimum (i.e., ones with ci > 0) may lead to extensive 
configurational entropy that might enable extremely slow 
dynamics, dll, [ill 

The simple key idea regarding "spatial chaos" is as fol- 
lows. It is well known that nonlinear dynamical systems 
may have solutions that exhibit chaos. This has been 
extensively applied in the time domain yet, formally, the 
differential equations governing the system may deter- 
mine not how the system evolves as a function of the time 
t but rather how fields change as a spatial coordinate (x) 
[replacing the time (t)]. Under such a simple swap of 
t X, we may observe spatial chaos as a function of the 
spatial coordinate x. In general, of course, more than one 
coordinate may be involved. The resultant spatial con- 
figurations may naturally correspond to amorphous sys- 
tems and realize models of structural glasses. A related 
correspondence in disordered systems has been found in 
random Potts systems wherein spin glass transitions co- 
incide with transitions from regular to chaotic phases in 
derived dynamical analogs. [TSj 

In the translationally invariant systems that form the 
focus of our study, an effective free energy of the form 

^ j dHS\x)-lf (50) 



is generally associated with single component (n = 1) 
systems of the form of Eqs. ( |4]). In Eq. ((50)) . [i rep- 
resents the deviation from the transition temperature in 
Ginzburg-Landau theories (or equivalently related to Eq. 

m)- 

Euler-Lagrange equations for the spins S(x) are ob- 
tained by extremizing the free energy of Eq. (|50p. 
These equations are, generally, nonlinear differential 
equations (as discussed in Appendix [B|) . As is well ap- 
preciated, however, nonlinear dynamical systems may 
exhibit chaotic behavior. In general, a dynamical sys- 
tem may, in the long time limit, either veer towards a 
fixed point, a limit cycle, or exhibit chaotic behavior. 
We should therefore expect to see such behavior in the 
spatial variables in systems which are governed by Euler- 
Lagrange equations with forms similar to nonlinear dy- 
namical systems. Upon formally replacing the temporal 
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{h)S(x) vs S'(x) 
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(c)5(x) vs S"(x) 
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(d)5(x) vs S"'{x) 
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(f)S"(x) vs S"'{x) 



S"(x) 

(g)5"(x) vs S<'<{x) 



FIG. 8: Classiness in system with as in Eq. (|49|l with 
Ci = 5, C2 = 4 and u = 1 and /i = 1 in Eq. (fSO]) . 



coordinate by a spatial coordinate, chaotic dynamics in 
the temporal regime map onto to a spatial amorphous 
(glassy) structure. 

In Fig. 8(a)[ we illustrate the spatial amorphous glass- 
like chaotic behavior that a one-dimensional rendition of 
the system of Eq. (gH]) exhibits. In Figs. |8(b)[j8(g)| we 
provide plots of the spatial derivatives of different order 
vs each other (and S(x) itself). 

Another example comes from the spatial analog of dy- 
namical systems with nonlinear "jerks" . It is well known 
that systems with nonlinear "jerks" often give rise to 
chaos [7^ "Jerk" here refers to the time derivative of a 
force, or, something which results in a change in the 
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(a)5(x) vs X 



S(x) 

{h)S'{x) vs S{x) 





S'(x) 

(c)5"(x) vs S'(a;) 



S"(x) 

(d)5(x) vs S"{x) 



FIG. 9: Example of aperiodic structure inspired by system 
with nonlinear jerks. Here J{S{x),S' {x),S" {x)) = -2S'{x) + 
{\S{x)\- 1) and initial conditions are 5(0) = -1, 5"(0) = -1, 
S"'(0) = 1 (chosen from Ref. jlil). 



acceleration of a body. Translating this idea from the 
temporal regime to the spatial regime, one can expect to 
obtain a aperiodic/glassy structure in a system for which 
the Euler Lagrange equation, Eq. (|B1[) may seem simple. 
For example, if we have the following, Euler Lagrange 
equation for a particular one-dimensional system, 



S"'{x)=J{S{x),S'{x),S"{x)), 



(51) 



with a non-linear function J{S{x), S'{x), S"{x)) then the 
system may have aperiodic structure. An example is de- 
picted in Fig. (0). 

We now discuss 0{n) systems and illustrate the exis- 
tence of periodic solutions (and absence of chaos) in a 
broad class of systems. 

The Euler-Lagrange equations for the system in Eq. 
([50)) [written longhand in Eqs. (|B1[ IB7|) ] become linear 
in case of "hard" spins, i.e., when the 0{n) condition is 
strictly enforced, i.e., u — >■ oo. In this limit, all configu- 
rations in the system can be described by a finite set of 
modulation wave- vectors (as was the case for the ground 
states in Sec. IV A|) . 

There are several ways to discern this result. First, 
it may be simply argued that since the Euler-Lagrange 
equations represent a finite set of coupled linear ordinary 
differential equations, chaotic solutions are not present. 
The configurations, therefore must be characterized by a 
finite number of modulation wave-vectors. 

A second approach is more quantitative. The idea used 
here is the same as the one used in Ref. [i^ . An identical 
construct can be applied to illustrate that spiral/poly- 



spiral states are the only possible states that satisfy the 
Euler-Lagrange equation if n > 1. With v being a func- 
tional of the lattice Laplacian of Eq. ([7|), the lattice ren- 
dition of the Euler-Lagrange equations in Fourier space 
reads 



D{Aj:)sik) = 0. 



(52) 



In what follows we consider what transpires when the 
Euler-Lagrange equations have real wave-vectors /C = 
{(fmjvas solutions. 



k — Qn 



0. 



(53) 



To obtain a bound on the number of wave-vectors that 
can be used to describe a general configuration satisfying 
the Euler-Lagrange equations, we consider general situa- 
tions wherein (i) 2{q,n±qm') 7^ Kec for any (f™, q,n' G /C, 
where krec represents a reciprocal lattice vector; and, (ii) 
qra ± qni' 7^ (fp ± <fp' for any qm.qra' ^ IC. Let a 
particular state be described as 



(54) 



where the vectors cim have n components for 0{n) sys- 
tems. As the states must have real components, the 
above equation must take the form. 



Soix!) = J2 



-tqm-x 



(55) 



In the above, a*^ denotes the vector whose components 
are complex conjugate those of the vector dm- In Eq. 
(j55p . wc do not count terms involving the wave- vectors 
(f,„ and —qm separately as such terms has been explicitly 
written in the sum. 

We next define the complex vectors {[/,„} and {V„i} 

as 



^iqm-x 



(56) 



The 0(n) normalization condition can then be expressed 
as, 

in 
m 



U*-Vm' + u*,-v„ 



(57) 



Solutions to Eq. ([57)) are spanned by the set of mutu- 
ally orthonormal basis vectors {Um} U {Kn}- As these 
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2Nq basis vectors are described by n-componcnts each, 
it follows that 

Nq < n/2. (58) 

Therefore, such states satisfying the Euler-Lagrange 
equations for an 0{n > 2) system can at most be charac- 
terized by n/2 pairs of wave- vectors. These states can be 
described by Ng spirals (or "poly-spirals" ) each of which 
is described in a different orthogonal plane. 
A few remarks are in order. 

• When u in Eq. ([50| is finite, i.e., in the soft spin 
regime, poly-spiral solutions could be present even 
though aperiodic solutions are also allowed. 

• Continuum limit: In the hard-spin limit, i.e., u — >■ 
oo in Eq. (jSOp , if the Fourier space Euler-Lagrange 
equation is satisfied by non-zero real wave-vectors, 
we have poly-spiral solutions as in the lattice case. 
When u is finite, aperiodic solutions may also be 
present. 

• If the Fourier space Euler-Lagrange equation docs 
not have any real wave-vector solution, poly-spiral 
states are not observed. 

In nonlinear dynamical systems, chaos is often ob- 
served via intermittent phases. As a tuning parameter 
A is varied, the system enters a phase in which it jumps 
between periodic and aperiodic phases until the length 
of the aperiodic phase diverges. This divergence is char- 
acterized by an exponent u = 1/2 similar to ours.[75| 



X. CONCLUSIONS 

Most of the work concerning properties of the correla- 
tion functions in diverse arenas, has to date focused on 
the correlation lengths and their behavior. In this work, 
we examined the oscillatory character of the correlation 
functions when they appear. 

We furthermore discussed when viable non-oscillatory 
spatially chaotic patterns may (or may not appear); in 
these, neither uniform nor oscillatory behavior is found. 
Our results are universal and may have many realiza- 
tions. Below, we provide a brief synopsis of our central 
results. 

1. We have shown the existence of a universal modu- 
lation length exponent = 1/2 [Eq. (|17p ]. Here 
the scaling could be function of some gen- 
eral parameter A such as temperature. This is ob- 
served in systems with analytic crossovers includ- 
ing the commensurate-incommensurate crossover 
in the ANNNI model. 

2. In certain situations the above exponent could take 
other rational values [Eq. [T5] . 



3. This result also applies to situations where a cor- 
relation length may lock in to a constant value as 
the parameter A is varied across a threshold value 
[as in Eq. (^1)) ]. 

4. We extended our result to include situations in 
which the crossover might take place at a branch 
point. In this case irrational exponents could also 
be present. In Eqs. p2[ I33|) . we provide univer- 
sal scaling relations for correlation and modulation 
lengths. 

5. We illustrate that discontinuous jumps in the mod- 
ulation/correlation lengths mandate a thermody- 
namic phase transition. 

6. We showed that in translationally invariant systems 
(with rotational and/or reflection symmetry), the 
total number of correlation and modulation lengths 
is generally conserved as the general parameter A 
is varied. 

7. Our results apply to both length scales as well as 
time scales. We further introduce the notion of a 
Joscphson time scale. 

8. We comment on the presence of aperiodic modu- 
lations/amorphous states in systems governed by 
nonlinear Euler-Lagrange equations. We illustrate 
that in a broad class of multi-component systems 
chaotic phases do not arise. Spiral/poly-spiral so- 
lutions appear instead. 

9. Our results have numerous applications. We dis- 
cussed several non-trivial consequences for classi- 
cal system in the text. For completeness, in Ap- 
pendix we discuss, rather simple applications 
of our results to non-interacting Fermi systems. 
We mention situations in which the Fermi surface 
changes topology, situations with divergent effec- 
tive masses and the metal-insulator transition. We 
further discuss applications to many other systems 
including Dirac systems and topological insulators. 
Aside from uniform and regular modulated peri- 
odic states of various strongly correlated electronic 
systems, there are numerous suggestions and 
indications of glassy (and spatially non-uniform or 
chaotic) behavior that naturally lead to high en- 
tropy in these systems, e.g., see, e.g., Refs. f72l.[76l- 
[roj . When spatial modulations are present in the 
ground states of rotationally invariant (and other) 
systems, they may lead to "holographic" -like en- 
tropy (as in large n renditions), |4l| . In future 
work, we will elaborate on non-trivial consequences 
of our results for interacting Fermi systems. 

Our general analysis regarding the expansion of the in- 
verse correlator as a function of k about points fc* 
and the myriad conclusions that we draw from it (includ- 
ing exponents) may, in some cases, be viewed as a formal 
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analog of Ginzburg-Landau method of expanding an ef- 
fective free energy T in an order parameter field (i.e., 
(5fc <-> and o T). 

Finally, we make a brief parenthetic remark concern- 
ing the "fractal dimension" in glasses and other systems. 
The notion of fractal dimensionality was recently applied 
in Ref. [13] based on a comparison between the atomic 
volume and the reciprocal of the dominant peak Kb, in 
the structure factor in metallic glasses. Specifically, the 



volume V 



with D f being the fractal dimension. 



This definition is very intuitive and such a relation be- 
tween volume and structure factor peaks is to be expected 
for a system of dimension Df if all natural scales in the 
parameter expand or contract with temperature (or other 
parameters) in unison. However, as we elaborated on at 
length, aside from global changes in the lattice constant, 
Kji can change non-trivially with temperature and other 
paramters in some regular lattice and other systems. For- 
mally, this may give rise to an effective non-trivial fractal 
dimension in various systems. 
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Appendix A: Fermi systems 

In this appendix, we discuss several examples of non- 
interacting fermionic systems where we observe a corre- 
lation or modulation length exponent. We will, in what 
follows, ignore spin degrees of freedom which lead to sim- 
ple degeneracy factors for the systems that we analyze. 
In non-interacting Fermi systems, the mode occupancies 
are given by the Fermi function. That is. 



(n(fc)) = {cHk)cik)) = 



g^)(e(fe)-M) 4-1' 



(Al) 



where c{k) and c^(fc) are the annihilation and creation 
operators at momentum k and /3 = l/(fcBT) with T the 
temperature. The correlation function associated with 
the amplitude for hopping from the origin to lattice site 
X is given by 



G(a-) = (Ct(G)C(x))=^(n(fe))^ 



-ik-: 



(A2) 



Thus far, in most explicit examples that we considered 
we discussed scaling with respect to a crossover temper- 
ature. In what follows, we will, on several occasions, 
further consider the scaling of correlation and modula- 
tion lengths with the chemical potential ^. We will use 



the letter v to represent exponents corresponding to scal- 
ing with respect to /j, and continue to use v to represent 
scaling with respect to the temperature T. 

The existence of modulated electronic phases is well 
known. [l-[ll [H, [H [H HI] In particular, the Fermi 
wave-vector dominated response of diverse modulated 
systems as evident in Lindhard functions, particular 
features of charge and spin density waves dominated 
by Fermi surface considerations in quasi- one dimen- 
sional and other systems have long been discussed 
and have numerous experimental realizations in diverse 
compounds. [8lll82j The exponents that we derived in this 
work appear for all electronic and other systems in which 
a crossover occurs in the form of the modulations seen in 
charge, spin, or other degrees of freedom. Our derived 
results concerning scaling apply to general interacting 
systems. To highlight essential physics as it pertains to 
the change of modulations in systems of practical impor- 
tance, we briefiy review and further discuss free electron 
systems. 



1. Zero temperature length scales — Scaling as a 
function of the chemical potential ^ 

We first consider a non-interacting fermionic system 
with a dispersion e{k). At zero temperature, the number 
of particles occupying the Fourier mode k is given by 



1 for e{k) < 1^1 
for e{k) > fi. 



(A3) 



All correlation functions as all other zero temperature 
thermodynamic properties, are determined by the Fermi 
surface geometry. We now consider the correlation func- 
tion of Eq. (jA2p . This correlation function will generally 
exhibit both correlation and modulation lengths. To ob- 
tain the modulation lengths along a chosen direction (the 
direction of the displacement x), a ray along that direc- 
tion may be drawn. The intercept of this ray with the 
Fermi surface provides the pertinent modulation wave- 
vectors. As we vary /i wc alter the density, p via 



(fc)<M 



d'^k 
(2^' 



(A4) 



gs being the spin degeneracy {gs = 2 for non-interacting 
spin-half particles such as electrons). As the Fermi sur- 
face topology is varied, the following effects may be ob- 
served. 

1. If two branches of the Fermi surface touch each 
other at p = po and are disjoint for all other 
values of ^, then a smooth crossover will appear 
from one set of modulation lengths to another with 
\Liy — Ldq\ oc l/i— /iol on both sides of the crossover. 
This crossover will be associated with an exponent 
vl — 1 characterizing the scaling of the modulation 
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|i for metal 



CO 




FIG. 10: Transition from a metal to a band insulator, 
figure is for illustration only. 



This 



lengths with deviations in the chemical potential. 
An example where a crossover of this kind is real- 
ized is the Cg = case of the schematic shown in 
Fig. [10] in which the crossover occurs at yu = fiQ. 
Other examples of this occur at half filling of the 
square lattice tight binding model and at three- 
quarters filling of the triangular lattice tight bind- 
ing model. These will be discussed later. 

2. If on the other hand, one branch of the Fermi sur- 
face vanishes as we go past /i = /io, the crossover 
is not so smooth and we get some rational frac- 
tion (usually vl = 1/2) as the crossover expo- 
nent: \Ld — Luq\ (X l/i — Hol'"'^ , on one side of the 
crossover. An example of this is shown in Fig. [TT] 
Here, 



\Ln-L 



D0\ 



2\n - fio\ 



2tt V |e"(27r/LDo)r 



(A5) 



where Ldq is the modulation length at the point 
where the ^ = /io line touches the e{k) curve, such 
that e'{2'!r/Luo) = The hopping correlation func- 
tion takes the form. 



G{x) 



2™d 



(cx)'^/2j^/2(ca:) 



(27r)d/2a;<i 



{2nY/ 



(A6) 



where /ig < /i < fiQ and a, b and c in Eq. (jA6[) 
(corresponding to modulation lengths of 27r/a, 27r/6 
and 27r/c) are the values of k for which e{k) = fi 
(as shown in Fig. [TT|) . 

At arbitrarily small but finite temperatures, the corre- 
lation function exhibits modulations of all possible wave- 



lengths. The prefactor multiplying a term with spa- 
tial modulations at wave-vector k is the exponential of 
{—\e{k) — An illustrative example is provided in 

Fig. 1121 Apart from the dominant zero temperature 
modulations, associated with the wave-vector ^2 in Fig. 
I12[ at finite temperature, there are additional contribu- 
tions from wave-vectors for which |e(A:) — fi\ is small rel- 
ative to fcsT. Near k2, we can assume e(fc) is linear 
such that e(k) w ii + (k ~ fc2)e'(fc2)- Similarly, near fci, 
e{k) — ^ w —A — {k — ki)^e"{ki)/2, where A = fi — fio (see 
Fig. [T^ . For large (3, both these contributions are highly 
localized around ^2 and ki respectively making the above 
approximations very good and the Fourier transforming 
integrals easy to evaluate {{n(k)) taking exponential and 
Gaussian forms). We have, 



G{x) 



ik2x)''^^Jd/2ik2x) 2(fc2x)^/2jrf/2_i(fc2x) 



(27r)d/2a;d 



{2TrY/2/3e'{k2)x'i- 



-/3A 



{kixf^^3d/2Mhx) 



(27r)— y^^^r 



(A7) 



where /3 — oo and A — )- 0, such that /?A — oo. 

Next, we will discuss scaling of the modulation length 
in with the chemical potential, /i in the familiar tight 
binding models on the square and triangular lattices at 
zero temperature. 



a. Tight binding model on the square lattice 

We consider a two-dimensional tight binding model of 
the square lattice. The dispersion in this model is given 
by 



(A8) 



£{k) = — 2t (cos kx + cos ky) 



The constant energy contours corresponding to Eq. 
are drawn in Fig. [131 As is clear from Fig. [131 there are 
certain directions (e.g., along the X-axis) along which 
there is no k for e{k) > 0. If we consider the same system 
at zero temperature, the following three crossovers are 
observed. 

(i) Half filling: The chemical potential ^ is zero at the 
half filling state. The Fermi surface is given by zkk^iky = 
TT. For small fi, we have, 



± kx i ky ^ TT -\ 

2t sm 



(A9) 



thus giving us an uninteresting modulation exponent, 
Vl = 1- 

(ii) Empty band: When fi = — 4i, none of the states 
are occupied. As we increase ^ hy a, tiny amount 6fi 
above this value, we observe a non-zero modulation wave- 
vector, k ~ y^Sfi/t, thus showing a modulation exponent 
VL = 1/2. 

(iii) Full inert bands: When /t — +At, all the states are 
occupied. As we lower /i by a tiny amount 5^, below this 
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(a) 



(b) 



FIG. 11: Example of a Fermi system where the modulation length exponent is 1/2. The gray region shows the filled states. 
When fj. > Ho, modulations corresponding to wave- vectors k = a and k — b cease to exist and we get an exponent of 1/2 at 
this crossover. Similarly, when h < fiQ, modulations corresponding to wave- vectors k — h and k = c die down. 




FIG. 12: The same Fermi system as in Fig. 1111 but now with 
a chemical potential fi = jj,o + A, slightly higher than fio. The 
temperature is small but finite. 



value, we observe a difference 6k of the modulation vector 
from ±ea;7r ± CyTT. We have, 6k = y^^Ji/t, thus showing 
a modulation exponent vl = 1/2 again. 



b. Tight binding model on the triangular lattice 

The analysis of the triangular lattice within the tight 
binding approximation, is very similar to the square lat- 
tice discussed above. The dispersion e(fc) is given by 



/I \ i~i 1 i kx ky'\/^ 

e[k) = —2tcoskx — 4i cos — cos — — . 



(AlO) 



We have exponents similar to the square lattice. 

(i) Three-quarters filling: The chemical potential fj, 



2t corresponds to the three-quarters filling state. If we 
concentrate on the {kx = tt, ky : — )- tt/VS} 
segment (same phenomenon is present at all the other 
segments of the quarter filling Fermi surface), we get, 



6kx 



6fi 



2 cos 



(All) 



where kx = + 6kx is obtained when fj. = 2t + 6fi. This 
leads to a modulation exponent oi vl — 1- The Fermi 




FIG. 13: Constant energy contours for two-dimensional tight 
binding model on the square lattice in Eq. (|A8|I . The red 
square corresponds to the particle hole symmetric contour 
where e{k) = 0. The contours inside it are for negative e{k) 
and those outside are for positive e(fc). 
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surfaces for chemical potentials /i close to three-quarters 
filling are schematically shown in Fig. [T31 

(ii) Empty band: When /i = — 6t, none of the states 
is occupied. As we increase ^ by a tiny amount 6fi 
above this value, we observe a non-zero modulation wave- 
vector, k = y^cJ/I/S, thus showing a modulation expo- 
nent vl ~ 1/2. 

(iii) Full inert bands: When fi ~ 3t, all of the states 
are occupied and close to this value the Fermi surface is 
composed of six small circles around k = a; cos(ri,7r/3) 4- 
y sin(n7r/3), n = {0, 1, 2, 3, 4, 5}. E ~ 3t — 6fi, we get, 
|5fc| = 2y^ (5/x/3, again giving us a modulation length ex- 
ponent, vl = 1/2. 



c. Metal-Insulator transition 

We discuss here the metal to band insulator transition 
at zero temperature. In a non-interacting system, this 
occurs when the Fermi energy is changed such that all 
occupied bands become completely full, as shown in Fig. 
1101 In the insulator, the Fermi energy lies in between 
two bands and thus the filled states are separated from 
the empty states by a finite energy gap. As the Fermi 
energy is tuned, the Fermi energy might touch one of the 
bands thereby rendering the system metallic. Close to 
this transition, the energy is quadratic in the momentum 
k, i.e., |fc| (X This implies that. 



|<5fc| cx \6^l\^/\ 



(A12) 



Following the scaling convention in Eq. ([T7|) . we adduce 
a similar exponent 



VL = 1/2 



(A13) 



that governs the scaling of the modulation lengths with 
the shift dfi of the chemical potential (instead of temper- 
ature variations). 



d. Dirac systems 

The low energy physics of graphene and Dirac systems 
is characterized by the existence of Dirac points in mo- 
mentum space where the density of states vanishes and 
the energy, e(fc) is proportional to the momentum k for 
small k. When we invoke and repeat our earlier analysis 
to these systems, we discern a trivial exponent 



\6k\ oc 



VDi 



= 1. 



(AM) 



This exponent may be contrasted with that derived from 

Eq. (imi . 



e. Topological Insulators - Multiple length scale exponents 
as a function of the chemical potential fj, 

The quintessential low energy physics of three- 
dimensional topological insulators can be gleaned from 
the following effective Hamiltonian [ssj in momentum 
space. 



iJ(fc) 



eo(fc)/4x4 + 
/M{k) A,k, 



Aik, 


V A^k^ 



^M{k) 
A2k+ 






M{k) 
-Aik, 



A2k. ^ 



-Aik^ 
-M{k) J 



A15) 



where eo(fc) = C + Dikl + D2k\, M{k) = M - Bikl 



B2k^ , with k± ^ kx + iky, k 



1,2 



y, n,j_ ^ y K- -r fc2 and Ai, A2, 

Bi, B2, C, Di and D2 constants for a given system. The 
energy bands are given by 



e(fc) = eo(fc) ± \/M{k)^ + Aikl + Aafc^. 



(A16) 



These bands are plotted in Figs. 15(a) and |15(b)] The 
finite gap between the two bands leads to an exponen- 
tially damped hopping amplitude, characterized by a fi- 
nite correlation length when the Fermi energy lies within 
this gap. These energy bands disperse quadratically for 
small k thus yielding 



\5k\ 



Vbuik = 1/2 



(A17) 



whenever the correlation length diverges and a insulator 
to metal transition takes place in the bulk, thus allowing 
long range hopping. The same exponent is also expected 
whenever the modulation length becomes constant as /i 
crosses some threshold value. 

The effective Hamiltonian for the surface states is given 

by 



H. 



surf 



Aafc^ 
A2k+ 



leading trivially to surface energies 

^surf{kx, ky) ~ ±A2kj_. 



(A18) 



(A19) 



Similar to the Dirac points in graphene (see Fig. 15(c) I, 
we trivially find an exponent of 



1. 



(A20) 



/. An example of a zero temperature Fermi system in which 
vl is not half or one 

Very large (or divergent) effective electronic masses 
TOe/ / can be found in heavy fermion systems (and at pu- 
tative quantum critical points). (s^. [Sq If the electronic 
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FIG. 14: Fermi surface for a triangular lattice with tight binding. The dashed lines are the Brillouin zone boundaries. This 
demonstrates a smooth crossover from one set of Fermi surface branches to another as fj, is changed across fi = 2. The points 
where the crossovers take place are (0, ±27r/v^), (±7r, ±-K/y/Z). The modulation length exponent for this crossover is vl = 1. 
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FIG. 15: Energy levels of Bi2Sez topological insulator. 



15(a) e(fc) versus k± at = 0; 15(b) e(fc) versus kz at k±_ = 0; 15(c) esurf{kx, ky) versus k± = {k^, ky) 



dispersion e{k) has a minimum at fco then a Taylor ex- 
pansion about that minimum trivially reads 

^Aiji{h - kQi){kj - kQj){ki - koi) + .(.A21) 

ijl 

When present, parity relative to fco or other considera- 
tions may limit this expansion to contain only even terms. 
As an example, we consider the dispersion 



e(fc) = ci - C2(fc^ 



(A22) 



where C2 > 0. The hopping correlation function of such a 
system has a term which exhibits modulations at wave- 
vector k ~ kQ at II = fi^, ~ ci. At higher values of the 
chemical potential, such a term ceases to exist. At lower 
values (fi = fi^, — Sfi) , this term breaks up into two terms 



whose modulation wave-vectors are different from fco by. 



fc — fco 



2fcoC2^^ 

VL = 1/4. 



(A23) 



2. Finite temperature length scales — Scaling as a 
function of temperature 

At finite temperatures, apart from the modulation 
lengths, there generally is a set of characteristic correla- 
tion lengths. From Eq. (jA2p . these are obtained by find- 
ing the poles (or other singularities) of the Fermi func- 
tion. Along some direction eo, the wave- vector fco = eofco 
is associated with a pole fco = ±27r/Lo ± */Co- At this 
wave- vector, 



e(fco) 



2n+ 1 



(A24) 
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where n is an integer. For a given /i, let us suppose that as 
we change the temperature, at T = Tq, we reach a saddle 
point of e(fc) in the complex plane of one of the Carte- 
sian components of k. Then, near this saddle point, the 
corresponding correlation and modulation lengths scale 
as, 

\ld-ldo\ cx iT-nr, 

le-Col « |r-ror% (A25) 

where = i^c = 1/2 in most cases (when the second 
derivative is not zero). 

Appendix B: Euler-Lagrange equations for scalar 
spin systems 

We elaborate on the Euler-Lagrange equations associ- 
ated with the free energy of Eq. ((50l) in Sec. |lXl These 
assume the form, 

J d''yV{x^y)S{y) + iiS{x) 

+u{S^{x) - \)S{x) = 0, (Bl) 

where V{x) ~ \V{x) + V{—x)]/2. For example, if we 
consider the finite ranged system for which, 

j d''yV{x~y)S{y) - aV'S{x) 

+ by^S{x) + (B2) 

then, we have, 

aW^S{x) + bW'^S{x) + ... + fiS{x) 

+ u{S^{x) - l)S{x) = 0. (B3) 

For lattice systems, the Euler Lagrange equation (jBip 
reads 

J2v{x~y)S{y) + ^lSix) 

V 

-f u{S^{x) - l)S{x) = 0. (B4) 
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In general, it may be convenient to express the linear 
terms in the above equation in terms of the lattice Lapla- 
cian A. We write 

D{^)S{x) = Y,y{^-y)S{y) + t^s{x), (B5) 

y 

D being some operator which is a function of the lattice 
Laplacian A. The real-space lattice Laplacian A, given 
by the Fourier transform of Eq. ([7]), acts on a general 
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